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' We present a model for the metal-insulator transition in 2D, observed in the recent years. Our 

I starting point consists of two ingredients only, which are ubiquitous in the experiments: Coulomb in- 

04 ■ teractions and weak disorder spin-orbit scattering (coming from the interfaces of the heterostructures 

^ ' in question). In a diagramatic approach, we predict the existence of a characteristic temperature 

^ , To — To{n,ujH), n being the density of carriers, and uh the Zeeman energy, below which these 

• systems become metallic. This is in very good agreement with experiments, and corroborates the 

fact that varying n and uh are equivalent ways into/out of the metallic regime. The conductivity, 
calculated as a function of temperature and loh in the metallic state, compares favorably to exper- 
iment. Moreover, we give an explicit expression for the conventional weak disorder contributions 
(Ih \ to the conductivity in the frame of our model. We comment on the nature of the transition, and 

Ch . calculate the specific heat of the system. 
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PACS numbers: 71.30.+h, 72.10.Bg, 72.10.-d, 72.15.Rn 



. - In the last years, a metal-insulator transition has been observed in 2D systems by Kravchenko et al.0 and others 
2 ' see e.g.l3. In high mobility heterostructures, for carrier densities n higher than a critical density ric ~ \Q^^cm~^ , a: 



■ at appropriately low temperatures T, of a few degrees Kelvin at most, a transition into a metallic state is observe' 
' Moreover, a sufficient increase of the Zeeman energy uiu of the carrier, through an exterrially applied magnetic field 

^ , H, takes the system back into the insulating statea. Hamilton et al.0 and Pudalov et al.cl have shown that beyond 

Q a maximum characteristic n'^ there is a second metal ta.jjisulator transition. Further, in the metallic regime, the 

O resistivity is usually fitted very well by the formula - e.g.Ero: 

^ ; Pe.p{T) = Po + Pi cxp [-{T^/Tf] , (f) 

Q.^ ; where po, pi ^nd depend on ri. The exponent k is in the range 0.5 - 1, and is material-dependent. 
^.f^ ' Hani et al.El have shown that for n < ric the insulating state is spacially inhomogeneous. Here, we will not attempt 
fS| , to provide a description of the insulating state. We only give a mechanism for the transition and a description of the 
0^ ■ metallic state. These are based on strong spin-dependent particle-hole correlations, which arise for the appropriate 
' range of n and ujh, in the frame of a Fermi liquid formulation. 

1 The experimental systems in question havelwo ubiquitous characteristics. The ratio of the Coulomb energy to the 
Fermi eneegy ep is typically in the range 4-400. The importance of interactions has been decisively demonstrated by 
Hani et al.Q and Dultz and JiangQ, who probed the compressibility in the metallic regime, and found it in agreement 
with a many-body interacting picture. On the other hand, weak disorder spin-orbit scattering is inadvertently present, 
^ , coming from the interfaces of the heterostructures in question. These two characteristics constitute our starting point. 
• We note that magnetic spin disorder would have the same effect as spin-orbit disorder, in the fr|ame of our model. 
! Strong support for the role of spin-dependent scattering is provided by the expts. of Vitkalov et al.la, which show that 
O ■ spin up and down mobilities remain comparable for all H. This is naturally interpreted in terms of strong spin up and 

■ down mixing through disorder scattering. Further, Hani et afl have demonstrated the existence of localized charged 
J> . islands in the metallic phase, which can be magnetic. If this is the case, we can assume that these magnetic islands 

sit on top of the aforementioned spin-orbit impurities. Then the scattering strength of the impurities is enhanced and 
the metallic behavior is more pronounced - c.f. e.g. eq. (6) for To{n,ujH) etc. 



' In the foregoing, we consider the coupling U between opposite spin carriers, and we denote by 

u = UNf , (2) 

the dimensionless coupling, Np being the density of states at the Fermi level. As usual, we work in the regime epT > 1 
{Ti — 1), ep being the Fermi energy and the total impurity scattering rate. As we have shown int2l, in the presence 
of weak disorder, which includes spin scattering, the ladder diagrams in the particle-hole channel give rise to the 
propagators A^{q,iLjm), j = 0, ±1, obeying the coupled Bethe-Salpeter equations 

+ UV^A^ + UV°A° , (3) 
A° = UV°A^ + UV-^A° . 
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The variables q,Wrm which stand for the momentum and Matsubara energy difference between particle and hole lines 
respectively, were supressed. are given by V^^ = P^-^^ V° = [V°'° - V^^°]/2, X'-'^™^ (q, w^) = NpiDq^ + 
jATg^ /Sy/lDq^ + jATg^ /i + \LOm \ — imjU)H}i with Tg^ the total spin scattering rate and D the diffusion constant. In 
Appendix A we present the derivation of these equations and their solution. l— , 
What turns out to be of interest here, is the "dynamic limit" Dq^ < cOm, where the solution of these eqs. igl3 

""^'^''""^^ A^Dq^+B^\u^\+C^ ' 

with 

Au = 12 - 20u + 15uV2 + 6f7|j , B„ = 4 - 6u + iu^ /2 + 2n]j , (5) 
C„ =r[ l-2u + 3mV4 + 172^(1 -iiV4) ] , K^o = Uru{l + , Kui = Ur{l-u) , 
LM^2Uu{2 + n]j) , Lui^W{l-u) , Muo^2Uu{2 + n]j) , Af„i = J7(4 - 3u) , 

where Q.h — and r = 4t^^/3 0. 

We consider first the case with luh = 0. For u = 2/3 — 2, Cu < 0. For u > 0.845, i?„ < as well. Then the ratio 
Cu/Bu is negative for u = 2/3 — 0.845, with C„ < and i?„ > 0. We interpret this as the onset of strong particle-hole 
correlations in the spin density channel. For sufficiently low temperatures 

r<T„(»,.„).J^.^^. (0) 

a resonance occurs - c.f. also eq. (^) , which we interpret as driving the second-order transition into the metallic 
regime - c.f. below. A finite ujh shrinks the range of u, which allows for this effect, until for high enough luh 
Cu becomes positive definite, and there is no transition into the metallic state. The dependence of-thfj-transition 
temperature To on u - or equivalently the carrier density - and loh are in accordaaee with experimenil3'Qll3i3. Namely, 
the lobe within which the conducting state exists in the ojh vs. density diagrarrfl'cl is easily reproduced by use of the 
ratio Cu/Bu- Of course, the value of u is determined through the bandstructure by the density n. 

We note that additional insertions in the Bethe-Salpeter equations (|^) for A-' , such as self-energy diagrams, do not 
influence the existence of this resonance. Based on the experimental evidence available to date, there is no magnetic 
or superconducting instability involved here. Thus-jthe location of the pole of on the imaginary energy axis is 
unphysical. In a different approach, Chamon et al.tj assume a paired state in the metallic phase. In order to make 
progress, we assume that in the metallic phase the usual Fermi liquid picture continues to apply, but C„ acquires a 
small imaginary factor 

Cu ^Cu + id , (7) 

which moves the pole of A^ at a distance off the imaginary axis, thus yielding well-deflned metallic contributions to 
the conductivity, d should be such that the contribution of the diagram B in fig. 1 is substantial for the calculation 
of the conductivity, as mentioned in Appendix B. Namely the quantity (c.f. eqs. ( p^ ) and (^6|) for Rj and eqs. ( [l0| ) 
and (|2|) for D{q,u;)) 

(2t)'^\R R' I 

= max{Gfl,(fc,e)GA(fc,e)}|i?,i?;i {coefiicient of ^ i?(<7, wt)} = \J^2ji ^ > (8) 
should not be too small. Also, d should satisfy 

d^±do\Cu\^T^ , (9) 

with < C < 10, 

< X ^ 1 and T^do{T) —> for T ^ 0, in order to yield physically meaningful contributions (the 
sign ± corresponds to the position of A^ with regards to the main particle-hole lines, as shown ia figs. 1-2). In fact, 
a finite d is consistent with the finite imaginary part of the dephasing rate for finite lo - c.f. e.g.cJ. 

To calculate the conductivity using the renormalized A^ , we proceed as follows. First, relevant processes must take 
into account the resonance of A^{q,u}). For this, diagrams, such as the ones in figs. 1-2, containing pairs of A^{q,Ld) 
's with identical q and lo axe needed in order to yield the enhancement factor Qjy of (j^) above, which then appears 
in eq. ([l|). The diffusonsEl 

D=^^' {q, u:ra) = 1^^^ 2?g2+j4r-V3+K|-zm,c.^ 
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in between the A^'s, both provide essential T-dependence and separate the 's from each other. We note that without 
something separating the 2 A^ 's, they would just reduce into a single A^ , and thus not offering this enhancement 
factor. The diagrams shown in figs. 1-2 yield the maximum coplribution around a single diffuson D{q,u!,n)- The 
reasoning here is similar to a related, but different, calculation inEJ. 

Summing over all possible spin combinations involving the different A^ 's yields for the total contribution of B 

B ^ g{ziFi + Z2F2} . (11) 
Here, g = mim {l^^T^Y ' ^1 = (^«l+A^«l^o)'' + (i^«o + Mu0^o)^ Z2 = 2(X„i +Af,at^o)^(i^«o + M„o'^o)^ , Fi = 

^I^l 27rT j ' ^2 - 2 in 1^ (27rT+r)2+c^;^ J and J i - T . 

Subsequently, we also take into account all the diagrams of the type of fig. 2, with one or two impurity lines passing 
completely outside the bubbles containing A^^s. Summing all diagrams from figs. 1 and 2 gives 

M = Ti - rr2 , (12) 



where T = Gnik, eF)GA{k, ep), Ti = ar^B, = B/{2TTeFT), with a = 1 - (2r)2 (glj + 3^ 

To calculate the conductivity, taking into account M to infinite order, we sum the infinite series, the n-th term of 
which contains n factors (blocks) M in series, sandwiched between the current vertices of a conductivity bubble : 

n, = '^ldkki{ry:{rM)" = -^^}. (13) 



m is the carrier mass. In this way we obtain a Af -dependent formula for the total conductivity in the metallic regime 

2NFe'^eFTT j y+ V- \ trr\ ^^^\ 

(T = (Jo + (TM = 5— — < , „ , „ ) Sa{T) . (14) 

3mS [ - - y_ J 

Here y± — (Fi ± S')/2, S = — 4r2. ccrnis the Drude term. Sa-{T) is a smooth analytic function, equal tOpl-for 
T <Ta and gradually vanishing for T > tJB. Eq. (|l|) is a gauge-invariant approximation for the conductivitycj'ta. 

The fit of eq. (|lj) to the activated exponential in T form for the resistivity of eq. (|l|) , usually fit to the experimental 
data, is shown in fig. 3. We take x = 1 " which gives an excellent match with eq. (|l|) in the "high T" limit, and makes 
the factor g T - independent - ^ = and do =const. We plot p(H ~ 0)/p*, where p — 1/a and — iraT / {"ne^ N ptp) . 
Sa{T) = 1 for all T here. The overall variation of p/p* within a given T range increases with u and g. Hence, 
an appropriate choice of|-t|hese parameters yields a resistivity which appears practically constant - c.f. e.g. the 
experimental data in ref.li3 - although it has the same shape as the curves of fig. 3. Overall, the fit of eq. (^^ is 
good. It can be seen that there is a discrepancy for low T . At the moment, we cannot determine the precise origin 
of this effect. 

In fig. 4 we plot the magnetoresistance corresponding to eq. ([l^), as a function of (H/H^). The field at which 
C„ = is -ff, = ^|1 — 2u + 3m^/4|/(1 — u'^ /i)/ {gcfiBr), gc the gyromagnetic ratio of the carriers and fiB the Bohr 
magneton. Again, Sa{T) — 1 for all T. The magnetoresistance for if scales like {H'^ /T'^), where 7 = 1.05 for the 
parameter set (1) of fig. 3. The typical range of values of 7 ~ 1 are in marked difference frpm — 2 expected from 
the conventional weak disorder magnetoresistance, but compare very favorably with expts.Ej. 

Turning to the conventional weak disorder conductivity contributions in the presence of finite M , it is easy to see 
that they can typically be written as 

2e^ f ^ T^W 

in a manner analogous to eq. ( |l3| ) - but note the square in the denominator of eq. (p^). W stands for any 
such diagrammatic contribution, with or without iateractions. E.g. the conventional weak localization correction 



corresponds to = — C, Cq,^^ being the Cooperonllj. From eqs. (13) for gm and (|T^) for aw, we obtain aw ^ (^m 
for W <C M . We believe this-|tO|be the explanation behind the negligible weak localization etc. contributions observed 
in a number of experimentstBEil. The same is true for the Hall coefficient as well, e.g.E2l. Note that, with increasing ii, 
M decreases, and the weak disorder aw is enhanced, as seen intiJ. On the othct-iiand, finite but small weak disorder 
contributions, which behave like -ln(T) in 2Dli3l3, have also been observed incBE^. In the frame of our model, it is 
reasonable to expect that at exponentially small T aw should dominate over ctm, and hence drive the system into 
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the insulating regime. In that case, we would not have a T = quantum metal-insulator phase transition. However, 
there are two caveats here. Kopietz has shown that the quasiparticle weight vanishes for T Oo, in an interacting 
2D Fermi gas with spinless disorder - the case with spin disorder remaining unresolved. Moreover, the well known 
Kohn-Luttinger instability to a superconducting phase may also interfere in the limit T — > 0. Hence the fate of the 
metallic state for T — > is an open question. 

Based on our model, we evaluate the specific heat SCy of the carriers, corresponding to the same processes yielding 
aM in cq. ([l^). The free energy is 

Fm = NflOotKIti C^~^/ )sf{T) , (16) 

with K = ^ g\z\F\ + 2:2-^2] and bCy = —T&^Fm I dT"^ ■ Here we took into account the diagrams of figs. 1 and 2. 
sf{T) has the same properties as Sa-{T) above. As before, this result should not be valid for T ^ 0. We note the 
pronounced H dependence of SCy ■ 

To summarize, we have shown that interactions and spin disorder can drive the second order metal-insulator 
transition observed in 2D systems. The transition arises from the onset of strong spin-density correlations, for a 
restricted range of the carrier density, the Zeeman energy and the temperature. We obtain good agreement with the 
experimentally determined dependence of the conductivity on these parameters. 

I have enjoyed discussions and/or correspondence with P. Kopietz, P.T. Coleridge, A. Hamilton, E. Kiritsis, S.V. 
Kravchenko, V. Senz, S. Vitkalov and D. Williams. I acknowledge the hospitality of V.S. Vassiliadis in Cambridge. 



APPENDIX A 

Equations (3), derived i n0, are standard Bethe-Salpeter equations in the spin-density channel, encompassing spin- 
scatering disorder and interactions. They are shown in fig. Al, and they form a complete system, containing all 
possible combinations of the interaction U and disorder lines in ladders, i.e. in the "no-crossing lines" approximation. 

Their solution is given by 

= C/(l - UV-^)/Do , A° = ~U^V^/Do , (17) 

where the determinant Do is 

Do = {l-UV-'^){l~UV^)~{UV°f . (18) 

To obtain the solution, we substitute I?' according to the expressions given in the text. We perform the polynomial 
expansions in q^, ujm, (jJh in both numerator and denominator of the expressions resulting from eqs. (0), after writing 
common denominators which are products of all relevant denominators coming from the I?"s. Then, takiiig the limit 
Dq^ < ujrn, yields the solutions in eq. (0). The solutions in the opposite limit Dq^ > uj,n can be found inEi 



APPENDIX B 



In this Appendix we calculate the basic diagrammatic block B of fig. 1 

B^T^ J2 G(fci,ei)G(fci -qi,ei -c^i)G(pi,ei)G(fc-(7i,e + r!-c^i) (19) 

{ei, 1^1,62, 1^2 }{fcl,Pl,<2l,fe2,P2, 92} 

G(fc2, e2)G(fc2 + 92, £2 + ^^2)G(p2, e2)G(fc + 92, e + ^^2) 
A\qi,iUi)A\k' -k + qi,e' -e + LUi) A^' {q2, iU2)A^' {k" - k + q^, e" - e + LU2) (^i - ^2, d - £2) • 

Here the indices /, mi depend on j and j' . All the energies are Matsubara energies, but we supress the usual indices 
m, n for the moment. The Green's function is G(fc, e„) = 1/ {«£„ ~ £fc + *sign(e„)/(2T)}. Fixing the external momenta 
and energies so that 

k' = k" = k , e' = e" = e , (20) 

the labels fc, e appear in both left and right external lines of the diagram. Thus the diagram can be inserted in a 
conductivity bubble, and calculate the conductivity accordingly. I.e. as a lowest order approximation - c.f. eq. ( p^ ) 
for the full formula - B can contribute to the conductivity through the formula 
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cji = — dk K{GR{k,eF)GA{k,ep)YB . (21) 

Eq. ( pO| ) implies directly that the two pairs of the interaction propagators A>{q,uj) - corresponding to the upper and 
lower bubbles - have identical arguments. This gives rise to the resonance factors referred to in the text - c.f. the 
discussion around eq. (^ and Rj below. Equivalently, taking identical arguments for each pair of A-' (q, a;)'s, yields 
eq. (ig|). 

In fact, the resonance in question selects identical arguments q, lo for each pair of A' (q, cl;)'s in diagrams of the type 
of B above, in which this is possible. In general, one can construct other diagrams, similar in structure to B, that 
is with pairs of A-'(q,w)'s connected to a common bubble, but in these instances it is impossible to obtain identical 
arguments q, w for the two (g, w) 's. These diagrams present no special interest, as they contain no resonance factors, 
and are ommitted. 

Also, eq. ( pO| ) implies that 

Pi^ki , _P2 = ^2 • (22) 

Moreover, we can ommit the gi , , 52 , "^2 dependence of the two Green's fimctions in the external part of the two 
bubbles. 

Taking all this into account we can write 

B^GR{k,e)GA{k,e)\RjRj,\Y,D''""i<i^^T) , (23) 

where 

i?j=r2 J2 iA^{q,LOrn)f G\k,en)G{k-q,en-u;m) , (24) 

0<w„, ,0<e„<w,„,fc,g 

and we take ujx as the smallest non-zero Matsubara energy difference 

iVT = 2ttT . (25) 

For the q summations we take Dq^ < t^^. Then, 

where the functions F1.2 appear after eq. (|ll|) and Fi corresponds to D°'°, while F2 to D^-^ - c.f. eq. (|l^ for D^-'"^ . 

APPENDIX C 

In this Appendix we evaluate the factor Rj. We see that 

R,=Tj2iA'iq,L^m)fF{q,L0ra) , (27) 

with 

F{q,iOm) = T J2 Glik,en)GAik^q,en-u;m) ■ (28) 

We can ommit the g-dependence in F{q,uJm), q < 1/VDt being much smaller than the relevant momenta of order 
kp in G(k,e), and then we evaluate 

J(fc,a;™)=T GUk,en)GAik, ) (29) 

0<c„<a)m 

1 f°° 

= 2~ deTiF{e-)[G%{k,e)GA{k,e~iuJm)-Gji{k,e + iujjn)GA{k,e)] . 
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npi^) is the Fermi-Dirac occupation function, and for low T we have 

r'" ( I 1 1 1 1 
2TTi I(k, ujm) ^ de \ } (30) 

J-oo - Cfc +it)'' e - tCi^rn - ~ it [c + ILUrn - (k + it) € - ek - it ) 

= ' ( ^ I + ' ln( (^^"^^)'+^^ U. 

LJm + 2f LeF + - efc + i< ep-f-k + it ujm + V (e^ - e^)^ + (w,„ + t)^ / J 

t = l/(2r). Subsequently 



FKO =iVF / defc /(fc,c^m) (31) 

1 1 r / ep + iujm + it^ 







- In H — 2(cjm + t arctan + 2t arctan — 

' ej, + {uj„i + tf 



In the limit 



of interest here, we obtain 



< i < ep- (32) 



^(^™) = • (33) 

With F{u!„i) in hand, we proceed to calculate 

YM = V(A^ (9,c^™))' = (AJ'(:r = ^)9^c^,„))' (34) 

1 )- 1 [Au {Kuj + Muj M)- Luj {Cu + Bu\uj„,\)]'^ ^ ^-".^ 



AttD I A'^ C„ + B„ I I + A„x 

where we restrict ourselves to the dynamic limit Dq^ < uJm- 



x=0 



1 



= —{AuLljX + 2Lu-j[Au{Kuj + Muj\uj^\) - L„j(C„ + Bu\uJm\)] \n{Cu + Bu\Lo„^\ + Aux)] , (35) 

is an irrelevant term, which does not yield a significant contribution compared to the pole of Y{ujm) - despite the 
logarithmic singularity - and will be ommited. 
Subsequently we evaluate 

n ^ ST v(. , , \ Fj^^m) [AujKuj + MujUJm) - L^jjCu + BuUJ.m)f , . 

~ 16in^d e%AuD 

where we approximated the uJm sum by the term most closely located to the pole Um = - cf. eq. (^) - as C„ < 
and Bu > 0, and we ommited the term in eq. ( ^4[ ) with x = cOm, as it yields a resonance at lower {To, ujo} than in eq. 
(i- 

* Current address. E-mail : kast@iesl.forth.gr 
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Fig. 1 

The basic (lovifest order) diagramatic block for M^. Dg^^ is the difFuson and A the renormalized propagator of eq. 

(4). All the appropriate spin structure diagrams are kept. 
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Fig. 2 

Variations on the diagram of fig. 1. Here the difFuson connects the bubble lines. One or two impurity scattering 

line(s) decorate the main upper and/or lower line(s). 




T(K) 

Fig. 3 

Fit of eq. ( |l4| ) - continuous lines - to the experimentally determined resistivity form of eq. - dashed lines. Here 
H = and = imr / {ire^ N pip) ■ X — ^ ^'^d C = 0. Both lower (1) and upper (2) curves have m=0.83, r = r^^ = 2 
K, and only differ by the ratio g'-^-'/s'"^'' = 2. The parameters of eq. (|l|) are (po/p*) = 0.392/0.519, 
{pi/pj) = 0.61/0.48, = 0.192/0.1 K and fc = 1/0.98, for curves (1) and (2) respectively. 
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Fig. 4 

Magnetoresistance from eq. (p^), for the parameter set (1) of fig. 3. iJ* is given in the text. The curves bottom to 
top correspond to T=0.3, 0.5, 1, 1.5 and 2 K. We note that p/p* = 1 implies M = 0. The fuU curves to H = H^, for 
T = 0.3 K are shown in the inset, the dashed Hne corresponding to the parameter set (2) of fig. 3. 
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Fig. Al 

The Bethe-Salpeter eqs. (3). The thin double line stands for U. The factors are given in the text. Note the 
explicit spin indices of the various components of the diagrams. 
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